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Abstract 

We present a method for the reconstruction of events containing hadronically decaying 
r leptons at collider experiments. This method relies on accurate knowledge of the r 
production vertex and precise measurement of its decay products. The method makes 
no assumptions about the r kinematics, and is insensitive to momentum loss along the 
beam direction. We demonstrate the method using e + e - —>• /i + /.4 _ r + r _ events fully 
simulated in the ILD detector. 
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1. Introduction 

Studies of final states including r leptons are of interest at current and future collider 
experiments; as an example, the dominant leptonic decay of the Higgs boson is to a r pair. 
A new generation of high energy particle colliders 0, i, i a is presently under study. 
A key scientific aim of these facilities is to measure the Higgs boson’s properties with 
great precision, important aspects of which involve measurements of the r final state. 
An example is the use of measurements of the r spin state to probe the CP nature of the 
Higgs boson. The detectors being designed for use at these accelerators will be equipped 
with vertex detectors providing unprecedented impact parameter resolution (see e.g. [5|), 
giving rise to intriguing possibilities in the reconstruction of relatively long-lived states 
such as the r lepton. 

We report on a method which uses a high-precision vertex detector together with 
other tracking and calorimetric detectors to fully reconstruct the kinematics of events 
containing hadronically decaying rs [i.e. decays in which only one v is produced) in an 
unbiased way. We outline previously used techniques for the kinematic reconstruction 
of single-i' r decays in section [2] In section [3] we define a new procedure which, in 
certain topologies, can fully reconstruct the r kinematics with significantly less stringent 
assumptions than previous approaches. This new method is then applied to e + e _ —> 
/i + /.i _ r + r _ events in sectionHJ and we conclude in section[5] 


2. Previous approaches to r pair reconstruction 

In the case of events containing a pair of r leptons each decaying to a single neutrino, 
the following method, which assumes knowledge of the rest-frame and invariant mass 
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of the r pair, but no knowledge about the r production vertex, is often used at lepton 
colliders ( e.g. jc|, 7j). The r-pair rest-frame can be assumed to be the centre-of-mass of the 
colliding beams (in the case of the e + e _ —>• t + t~ process), or the frame recoiling against 
particles produced in conjunction with the r pair, as in the case e + e _ —>• p + p,~T + T~. 
The r decay products are then boosted into the assumed r pair rest frame, in which the 
energy of the rs is defined by the assumed invariant mass of the r pair. The r mass then 
constrains the r momentum to be at a fixed angle to the momentum of its hadronic decay 
products, defining a cone around the hadronic momentum. The two cones in an event, 
one per r, have either 0, 1, or 2 intersections, corresponding to the possible solutions for 
the r momenta. 

At hadron colliders, the unknown net momentum along the beam direction results 
in less available information to constrain the event kinematics. The invariant mass of 
r pairs can be partially estimated using the invariant masses of visible decay products 
and the missing transverse energy, or by applying the approximation that the v from r 
decay is collinear with the visible r decay products [gj. Another approach is to combined 
the measured momenta of visible r decay products with constraints on the r mass and 
global event transverse momentum balance, resulting in an under-constrained system. 
The likelihood of the various r decay topologies allowed by the constraints can then be 
usedjo choose a best solution, or alternatively to associate a weight to each solution 
e.g. 

If the t production vertex is precisely known, the use of the impact parameters of 
the charged r daughters (“prongs”) brings additional information. The knowledge of the 
production vertex can come from the reconstruction of particles recoiling against the rs, 
or from a priori knowledge of the interaction point, if the size of the interaction region 
is sufficiently smaller than the impact parameters of the r decay products. The use of 
the impact parameter vectors of charged r daughters, without full reconstruction of r 
decay kinematics, in the analysis of Higgs boson CP properties has been described in e.g. 
2l [13j, while their use in fully reconstructing di-r systems of known invariant mass and 


momentum have been demonstrated in 1 


icting di 
2 15). 


3. Method 

In this section we present methods that can, under certain conditions, fully recon¬ 
struct a t without assuming that it belongs to a t pair of particular invariant mass or 
centre-of-mass frame. In section IT21 we consider the reconstruction of hadronic single 
prong final states, in which a single charged hadron is produced with zero or more neutral 
hadrons and a single neutrino. Such decays account for 49.5% of r decays. Multiprong 
hadronic r decays, in which three or more charged hadrons, zero or more neutral hadrons, 
and a single neutrino are produced, account for 15.3% of r decays, are discussed in section 
13.31 Leptonic decays of the r (35.2%) provide significantly less measurable information 
about its decay kinematics due to the production of two neutrinos, and are not further 
considered in this paper. 

The method relies on precise knowledge of the r production vertex and the charged 
prong trajectories, and on the reconstruction of any neutral hadrons produced in the 
decay. Constraints on the invariant mass and lifetime of each r, and on the overall 
transverse momentum in the event, are then used to determine the r momenta. 
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3.1. Tau production vertex 

The uncertainty on the r production position should be small in comparison to the 
decay length of the r and the typical impact parameters of its decay products. In final 
states in which the rs are produced together with more than one prompt charged particle, 
the production vertex can be directly reconstructed on an event-by-event basis using the 
tracks of these particles (e.g. the us in the process e + e _ —>• p + p~t + t~). The proposed 
linear electron-positron colliders [1, |2| have rather small interaction regions, which may 
be used as an additional constraint on the interaction point, although this is not done 
for the results presented in this paper. 

3.2. Single prong r decays 

3.2.1. Tau decay plane 

In the case of single prong hadronic r decays, the trajectory of the charged prong 
(helical in the usual case of a uniform magnetic field) can be used to define a plane 
(hereafter called the “track plane”) which contains two vectors: d, the vector between 
the reconstructed interaction point (IP, assumed to be the r production vertex) and the 
point on the trajectory closest to the IP (point of closest approach PCA); and p, the 
tangent to the trajectory at the PCA. In the case of linear trajectories of the r and of 
the prong between the PCA and the r decay vertex, the r momentum, and therefore 
also the sum of the momenta of the other decay products of the r (neutrinos and neutral 
hadrons), are constrained to lie within this plane. 

The difference between the reconstructed track plane and the true r decay plane 
(defined by the r and prong momenta) depends on the decay length of the t, the accuracy 
with which the IP position is known, the precision of the charged prong trajectory, and 
the extent to which the linear approximation of the r and prong trajectories near the IP 
is valicQ. 

3.2.2. Parameterisation of neutrino momentum 

In this section we describe the parameterisation of the unmeasured neutrino momen¬ 
tum q, based on the measured prong trajectory and neutral hadron momentum. 

Neutral particles are measured as clusters in the calorimeters, or as identified conver¬ 
sions of photons into e + e _ pairs within the tracker volume. The momentum to be asso¬ 
ciated to calorimeter clusters can be estimated by assigning the energy of the calorimeter 
cluster, a mass hypothesis (e.g. zero in the case of photon-like clusters, the Kp mass 
for hadronic clusters), and the direction of a straight line connecting the IP and the 
energy-weighted mean position of the calorimeter clustefl 

The three-momentum k of the neutral hadronic system can be decomposed into com¬ 
ponents perpendicular to and within the track plane: and k|| respectively. Since the 

r momentum lies within the track plane, the hadronic momentum perpendicular to the 


1 The error in this linear approximation scales as the ratio of the decay length of the r to the radius 
of curvature of the prong: for a prong with px = 10 GeV/c produced by a 50 GeV r of average lifetime, 
in a field of 3.5 T, this ratio is < 10 —3 . An iterative approach, in which a first iteration uses the helix 
parameters at the PCA to the IP, while later iterations use the helix parameters at the calculated r 
decay position, should reduce any sensitivity to the prong’s curvature. 

2 Alternative definitions are possible: for example the line connecting a first estimate of the r decay 
position to the identified start of the calorimetric shower. 
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Figure 1: Parameterisation of the neutrino momentum in the track plane (qy) in terms of the vectors 
h || and f, and angle -0. The track plane is defined by the vectors p and d. 


track plane must be balanced by the neutrino, so the perpendicular component of the 
neutrino momentum qj^ = — kj_- 

The component of the neutrino momentum within the track plane can completely 
generally be parameterised as 

q|| = Q • (cosi/j ■ h|| + sinip • f), El (1) 

where Q is the unknown magnitude of the in-plane component of the neutrino momentum, 
h|| is the component of the total hadronic momentum (h = p + k) in the track plane, and 
the unit vector f = f/|f|, where f = hy X (d X hy), is within the plane and perpendicular 
to h||. 

Four-vectors p , fc, and q corresponding to the three-momenta p, k, and q, can be 
defined by means of appropriate invariant mass assumptions. The invariant mass of the 
sum of four-vectors p, fc, and q must be equal to the r lepton mass m T . This constraint 
allows us to write an equation involving Q and qy (which in turn depends on i/i): 

ml = (,q + h ) 2 

= (^/Q 2 + ki + E h ) 2 - (Qq„ - k ± + h) 2 (2) 

where Eh is the energy of the hadronic system. This can used to solve for Q: 

Q = 77 -(- & ± Vb 2 - 4ac), (3) 


3 We define x to be a unit vector parallel to x. 
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where 

a = B 2 — C 2 , b = 2AB, c = A 2 - C 2 k ± 2 , (4) 

and 

A = m 2 - m 2 h - 2h • kj., B = 2q|| • (h - kj_), C = 2E h , (5) 

where rrih is the invariant mass of the hadronic system. 

There are in general two solutions of Q for each choice of ip, which are complex in 
unphysical regions. Such complex solutions are rejected. In the case of two real solutions, 
we choose to denote the one with a higher energy neutrino in the laboratory frame as the 
“first” solution, and the other as the “second”. The first solution at ip = a corresponds 
to the second solution at ip = 7r + a, and vice versa. Since the angle between the neutrino 
and hadrons is typically small in the laboratory (due to the large boost of the r), it is 
convenient to consider the solutions separately, each in the range —tt/2 < ip < 7r/2. 

For each real solution of Q, the corresponding r momentum can be calculated. In 
conjunction with the prong trajectory, this allows the decay length and proper decay 
time of the r to be calculated. In the simple case of k = 0 [i.e. no neutral hadrons in the 
r decay), one solution for Q corresponds to a negative r decay length (the intersection of 
the r trajectory with that of the prong is on the “wrong” side of the IP), and can therefore 
be discarded. More generally, a likelihood A that the extracted lifetime is consistent with 
that expected of the r can be expressed in terms of the measured decay length of the r 
in the laboratory L, its Lorentz boost (/3, 7 ), and the mean lifetime of the r (~ 87/im/c): 
A = exp (—L/(/3 • 7 • 87 pm)) if L > 0, and A = 0 otherwise]! 

3.2.3. Choice of ip 

To determine the value of ip, additional information is required. One possible ap¬ 
proach would be to follow a statistical approach, applying a weight to each possible 
solution based on the likelihood that its decay time is consistent with that expected of 
the well-known r mean lifetime, and/or that the reconstructed r decay kinematics follow 
the expected distributions. If a hypothesis is made as to the energy of the r, or on the 
invariant mass of a pair of rs, this can also help choose appropriate ip solutions. 

An alternative approach, followed in this paper, is to consider the environment in 
which the r has been produced. We consider the class of events in which one or more 
single-^ decaying r leptons have been produced in conjunction with zero or more well- 
measured particles, together with zero or more particles escaping along the beam-line 
(e.g. ISR photons). Examples of such processes at an electron-positron collider are two- 
fermion production e + e _ —» t + t~ , Higgs-strahlung e + e“ —> HZ —>• (r + T~)(p + p~), and 
its major irreducible background e + e _ —> ZZ —> (r + r _ )(// + /r _ ). In such events, the 
overall pt of the rs and other visible particles is balanced, while there may be non-zero 
net momentum along the beam-line due to ISR. A natural way in which to determine the 
ip angles in an event (one per r decay) is to choose that combination which minimises 
the magnitude of the missing transverse momentum (p™ lss ) in the event. This method 
has the advantage of making no assumptions about r kinematics or about undetected 
particles escaping along the beam-line - also indispensable for application at hadron 


4 A more sophisticated treatment would take measurement uncertainties into account, therefore al¬ 
lowing small negative decay length solutions. 
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collider experiments, but does require that no additional neutrinos are produced in the 
same event. In the case of several solutions for which p™ lss is consistent with zero, the 
reconstructed r lifetimes can be used to determine the most likely one, or each solution 
could be weighted by a likelihood based on its lifetime and/or r decay kinematics. 


3.3. Multiprong r decays 

In the case of multiprong r decays, the prongs can be fitted to a common vertex, 
whose position V corresponds to the point at which the r lepton decayed. Together 
with a well-known production position P, this gives an estimate of the direction of the r 
momentum, in the approximation that the r trajectory is lineaiJl. The r momentum p r 
can be written as p T = tr, where t is the magnitude of the r momentum, and r = V — P. 

The neutrino momentum q can then be written as q = p T h, where h is the 
measured momentum of the visible hadronic system. Defining appropriate four-vectors q 
and h for the neutrino and hadronic system respectively, and requiring that the invariant 
mass of the decay products is equal to the r mass m T , gives 


ml = (q + h ) 2 

= (|q| + Eh ) 2 — (q + h) 2 
= (\tr - h| + E h ) 2 - t 2 , 

where Eh is the energy of the hadronic system. This can be solved for f, giving 

t = -^-(6 ± \/b 2 — 4ac), 

where 


a = El - (r • h) 2 , b = (r • h)(m* + m 2 h ), c 


( TO r + 


m 


2 \2 
h) 


+ mzEt 


( 6 ) 

( 7 ) 

( 8 ) 


There are in general two possible solutions for t which satisfy the m T constraint, however 
they are not guaranteed to be real. The full tau decay kinematics can be calculated 
for each real solution, including the decay length and proper decay time of the r. In 
appropriate event topologies, the p™ ss of the event can be used to choose between these 
two solutions, otherwise the previously introduced lifetime likelihood could be used. 

Since the r leptons produced at high energy colliders are typically highly boosted, 
the opening angle of the multiprong jet is usually small. As a consequence, the precision 
with which the r decay vertex is reconstructed is significantly worse along direction of 
the r jet than along the perpendicular directions. In the example considered later in this 
paper, in which rs are produced in Higgs boson decays, the length of the major axis of 
the vertex position error ellipsoid is typically of order 100 pm, while the other two axes 
have lengths of order 2 pm. 

In this case, it turns out to be better to consider the major axis of the vertex error 
ellipsoid as a single trajectory, and to analyse the event using the procedure developed for 
single prong decays in the previous section. In this way, the r decay plane is defined by 


5 We note that the r momentum direction can be directly measured for any decay mode if the r 
decay length is sufficiently large for it to produce hits in the vertex detector. 
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the IP and the centre and major axis of the vertex ellipsoid, and the exact decay position 
along the major axis direction is left as a free parameter. The distance between the fitted 
and reconstructed vertex positions, normalised by the uncertainty on the reconstructed 
vertex position, can be used to define a vertex likelihood, which can be used when 
choosing between several possible solutions. 

4. Example application 

We apply this method to e + e~ —>• fi + H events, in which the Higgs boson decays 
to a pair of rs. Such events are usually selected by considering the mass recoiling against 
the di-/i system, which has a peak at the Higgs boson mass. However the distribution 
of this recoil mass has a long tail due to beamstrahlung and ISR, particularly at higher 
centre-of-mass energies well above threshold. If the invariant mass of the rr system can 
be directly reconstructed, a rather cleaner selection of such events should be possible. 
The full reconstruction of the r decay kinematics and the rr centre-of-mass frame also 
allows the best use of the spin information of the rs: the properties of the r decays can be 
used to define their polarimeter vectors , and the correlations between the r polarimeters 
can be used to measure the CP-nature of the Higgs boson (e.g. My 

Events were generated at a centre-of-mass energy of 250 GeV using WHIZARD v2.2.2jl6j], 
assuming a Higgs boson mass of 125 GeV, and including effects due to initial state radi¬ 
ation and beamstrahlung (by CiRCEl). The rs were decayed by tauola++ vl.l.4[l7j. 
Three scenarios were considered: both rs decaying to n^is- both decaying to 7r ± 7r°^ r ; 
and both decaying to (ai —> 37t ± );a The simplest hadronic r decay mode is ir^is, but 
accounts for only 11.5% of r decays, while the pis mode has the largest branching ratio 
of 26.0%@ 

The resulting events we passed through the GEANT4-based MOKKA simulation of the 
ILD_ol_v05 detector model Q. This model consists of a vertex detector with three 
double-layers of silicon pixel detectors, a silicon strip-based inner tracker and forward 
tracking disks, a large time projection chamber within a silicon tracking envelope, fol¬ 
lowed by highly granular calorimeters: a silicon-tungsten ECAL and scintillator-iron 
HCAL. These are placed within a solenoid producing a 3.5T magnetic field, surrounded 
by an instrumented iron flux return yoke. The simulated energy deposits within the ac¬ 
tive detectors were passed through the standard ILD digitisation procedures to simulate 
detector signals. 

4-1. Event reconstruction 

The standard ILD reconstruction software was used to reconstruct charged particle 
tracks and calorimeter clusters. For the results shown in this paper, the reconstructed 
tracks were associated to primary muons and charged pions from r decays on the basis of 
matching to the simulated particle directions. In a full analysis, muon and pion identifi¬ 
cation should be performed. This is rather simple task in detectors with highly granular 
calorimeters such as ILD, which should not present any great difficulties. Distributions of 
the impact parameter in the plane perpendicular to the beam do, and of its uncertainty, 


6 The vast majority of r —> 7% ir°u T decays proceed via the p f , so we use pv as shorthand for tt 7T° is T 
in this paper. 
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are shown in fig. [2] for both prompt p, tracks and 7r tracks from r decay. The two p 
tracks were fitted to a common vertex to measure the event IP using the LCFIVERTEx[18I| 
package, achieving a typical precision on the IP position of better than 3pm in all three 
dimensions. 

Clusters identified in the electromagnetic calorimeter by the GARLIC [l9j] and PANDORAPFA 
algorithms were treated as photon candidates. Each cluster was assigned a momentum 
with magnitude equal to the reconstructed cluster energy and direction parallel to a line 
joining the nominal IP and the energy-weighted mean position of the cluster. Clusters 
were associated to 7r°s and their parent rs based on matching the cluster momentum di¬ 
rection to the simulated photon directions. Events were rejected if not all the simulated 
p, 7 r and photons could be matched to reconstructed particles. Such cases are typically 
due to imperfect reconstruction algorithms or the interactions of particles in the tracker 
volume. 

In order to improve the effective photon energy resolution, clusters associated to a 
7T° were subjected to a constrained kinematic fit imposing the 7 r° mass by varying their 
energies. No special effort was made to correct invariant mass biases occurring when 
reconstructing two overlapping clusters from 7r° decay. Figure [2] shows the precision with 
which the visible invariant mass is reconstructed in r —> pu decays. The accuracy with 
which the reconstructed track plane approximates the true r decay plane is demonstrated 
in the bottom right plot of fig. [H in which the “track plane error” is defined as the angle 
that the true r momentum makes to the reconstructed track plane. 

The reconstructed tracks from multiprong r decays were fitted to a common vertex, 
again using the lcfivertex package. In three-prong decays via the a\, the ellipsoid 
describing the vertex position uncertainty typically has a major axis length of 100/im, 
and the other two axes around 2pm. 

4-2. Dependence o/p™ lss and the lifetime likelihood on 

Figure [3] shows the dependence of the event p™ ss (defined as the magnitude of the 
component of P = JT p; transverse to the beam-line, where the index i runs over 
the momenta p; of the rs and other reconstructed particles in the event) and lifetime 
likelihood on the two r decay angles ijr^ in two events, one in which both r —> the 
other in which both r —> pv. The distributions are shown for the case in which the first 
Q solution is used for each r (which is where the best solution was found in these two 
events). Four minima at which p™ lss ~ 0 are seen, one in each quadrant. The lifetime 
likelihood (shown for the same Q solution), together with the p™ lss at the minimum, 
allows the best minimum to be chosen. 

The MINUIT package in the ROOt[21| analysis framework was used to find the minima 
of the p™ lss distributions separately in each of the four quadrants bounded by 0 ± = 

0, ±7t/ 2, in each of the four Q-solution combinations. If no minimum was found, or 
either r had a negative reconstructed decay length at the minimum, the quadrant was 
rejected. Of the remaining possible solutions, the one with the smallest value of p™ ss 
was chosen. 

4-3. Results 

4-3.1. Single prong decays 

Of events in which reconstructed objects could be unambiguously matched to sim¬ 
ulated particles, no good solution was found in 0.4% (2.4%) of events in the n ± i'(pi;) 





Figure 2: Properties of the considered e“*~e — —> H events. Top: the measured 2-d impact 

parameter do [left] and its uncertainty [right], for prompt [i tracks, and ir tracks from r —»• pv decay. 
Bottom left: the difference between the reconstructed and true visible r mass in t pis decays. Bottom 
right: the angle between the reconstructed track plane and the true direction of the r momentum. 
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Figure 3: Example of two e~*~e —>• fi^p, t~*~t events in which both r —> [top], or both r —>■ 

[bottom]. Contours of p[left] and lifetime likelihood [right] as a function of the two angles are 
shown. Regions with darker contours (lower p'^' LSS : higher lifetime likelihood) are preferred. In the upper 
event no positive decay length solutions were found for < 0, while in the lower event, no physical 
solutions were found for 0-2. 
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channel. This inefficiency is due to cases in which the reconstructed visible mass is larger 
than the r mass, or in which no solution with positive decay length is found. In fig. [4] 
we show the distribution of the chosen solutions’ p™ lss , for events in which both rs were 
forced to decay to either or pv. A rather good minimum (p™ ss < 1 MeV) is found 
in the large majority of cases, but a fraction of events have a larger p™ lss . These events 
with large p™ lss are due to mis-reconstruction of the event, and make up a larger fraction 
of pv than 7r ± i' decays due to the relatively worse precision of the calorimeters compared 
to the tracker. Events in which the reconstructed p™ lss is larger than around 0.5 GeV/c 
show a significantly wider peak in the r — r mass distribution, as shown in fig. [I] The 
mass distribution in the case of n±v ( pv ) decays has a central peak with a Gaussian 
width er ~ 0.6 GeV (1.1 GeV), and non-Gaussian tails to both higher and lower values. 
Of 7T^~v [pv) events with p™ ss < 0.5 GeV/c, 74% (67%) lie within 3cr, and 95% (89%) 
within 10 GeV, of the peak position. 

The same figure also shows the reconstructed invariant mass of the two r system in 
events with small or large ISR/beamstrahlung energy (demonstrating that the method 
is independent of the total centre-of-mass energy and boosts along the beam-line), and 
on the smaller of the true r decay lengths in the laboratory (showing that events with 
longer decay lengths tend to be better reconstructed). Figure [5] shows the r pair mass in 
the e + e _ —> p, + p,~H event sample, and a second e + e^ —> p + p~r + T~ sample in which 
the Higgs boson contribution was effectively turned off. Clear mass peaks due to the 
Higgs and Z° bosons can be seen, the widths of which allow a rather clean separation 
between these processes in both the considered r decay modes. 

4-3.2. Multiprong decays 

Multiprong r decays, studied using the r — > (ai — > 3tt ± )v channel, are reconstructed 
using two methods. In the “vertex” method, the reconstructed vertex position is directly 
used to constrain the r direction, and the r mass constraint used to find two possible 
solutions for the neutrino momentum. No real solution is found in around a quarter of 
t decays. The event p™ lss is used to choose the best solution. 

In the “decay plane” method, the central position and major axis of the vertex posi¬ 
tion error ellipsoid are used to define the t decay plane, which is then used in the same 
way as for single prong decays. Candidate r solutions are required to have fitted and 
reconstructed r decay vertices in the same hemisphere. The event candidate solution 
with smallest p™ ss is selected. When using this method, no good candidate solution was 
identified in around 4.5% of events. Figure [6] shows a comparison of these two multi- 
prong methods. As well as having a much higher efficiency for identifying a solution, 
the “decay plane” approach results in significantly better mass resolution. The di-r in¬ 
variant mass distribution for multiprong ai decays obtained by using the “decay plane” 
method are compared to other decay channels in fig. [5j showing that multiprong decays 
are reconstructed with a similar precision as r —> n^v decays. 

In fig [7] we compare r pair mass reconstruction using the methods developed in this 
paper with alternative approaches. The simplest is to use the invariant mass of only 
the visible decay products, therefore ignoring the neutrino contribution. This clearly 
underestimates the invariant mass, and results in a rather wide distribution. The collinear 
approximation assumes that the neutrinos are parallel to the visible tau decay products, 
and further requires the event’s pt to be balanced, allowing the neutrino energies to be 
estimated. The mass distribution resulting from the collinear approximation is peaked 
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Figure 4: Properties of the considered e + e —>■ H events after full r reconstruction. The value 

of the event at the chosen solution [top left]. The r pair mass distributions in different ranges of 

^m(iss) right], the total ISR/beamstrahlung energy Ejgp^ [bottom left], and L m i n , the smaller of 
the two true r decay lengths [bottom right]. Except in the top left plot, the relative normalisations of 
the various contributions are as in the event samples. 
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Figure 5: The reconstructed mass of the r pair in e"*"e —»• r“*"r events in which the r pair was 

produced via a Higgs (left) or Z 0 / 7 * (right). 



Figure 6 : The reconstructed mass of the r pair in e+e - —»• //T/i - t + t — events in which both rs decay 
into three charged pions. We compare the results of two r reconstruction methods: directly using the 
reconstructed vertex position, and using the position and major axis of the vertex ellipsoid to define the 
t decay plane. 
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Figure 7: Comparison of the reconstructed r pair invariant mass using different methods: the impact 
parameter-based method developed in this paper, the collinear approximation, and only the visible decay 
products. 


at the correct value, however it is significantly wider than that achieved by the methods 
developed in this paper. The difference in resolution is particularly significant in the n ± iy 
decay channel, where the impact parameter method gives a very sharp distribution. 

5. Conclusion 

An approach to the reconstruction of hadronic rs has been presented, which can be 
used in events in which the r production vertex is well known, and the trajectory of its 
charged decay products is precisely measured. This method works only in events in which 
no undetected particles with significant pt are produced together with the rs. In contrast 
to other methods, no assumptions are made on the centre-of-mass frame or the invariant 
mass of a r-pair system, nor on the r energies. The method is insensitive to momentum 
balance along the beam-line, so this method should work on e.g. e + e~ —> p, + p~ H events 
produced both near, and well above, threshold. 

The example analysed in this paper, in which a r pair recoils against a p pair, is 
the case in which the event p™ lss is most precisely measured. The Higgs boson will 
more commonly recoil against a hadronic system. These events are more difficult to 
reconstruct, and will have a less precisely measured p™ ss . The use of a constrained 
kinematic fit, taking into account the resolution with which the various quantities are 
measured, should lead to improved results. 

This reconstruction technique can be used at lepton colliders, as demonstrated in 
this paper, provided the precision of the tracking detectors is sufficient, particularly in 
the estimation of charged particle trajectories near the interaction point. To be used 
at hadron collider experiments, the r decay products must be sufficiently well identified 
and reconstructed, and all particles produced in the same interaction unambiguously 
identified in order to accurately estimate the p™ ss associated to the interaction. This 
latter point suggests that cases in which rs recoil against a relatively simple system, e.g. 
/i — p as considered in this paper or a single high p T jet, should be the most promising 
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in the richer environment of hadronic collisions at high luminosity, in which multiple 
interactions typically occur in each event. The method is insensitive to the unknown net 
momentum along the beam direction inherent in hadronic collisions. 

If applied to Higgs boson decays into r leptons, the use of this method will allow clean 
separation of signal events from irreducible backgrounds from e.g. Z decays, and allow 
precise reconstruction of the r spin information, providing a sensitive probe to measure 
the CP nature of the Higgs boson. 
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